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In the past few years, the bone field has witnessed great advances in genome-wide as- 
sociation studies (GWASs) of osteoporosis, with a number of promising genes identified. 
In particular, meta-analysis of GWASs, aimed at increasing the power of studies by com- 
bining the results from different study populations, have led to the identification of nov- 
el associations that would not otherwise have been identified in individual GWASs. Re- 
cently, the first whole genome sequencing study for osteoporosis and fractures was 
published, reporting a novel rare nonsense mutation. This review summarizes the im- 
portant and representative findings published by December 201 3. Comments are made 
on the notable findings and representative studies for their potential influence and im- 
plications on our present understanding of the genetics of osteoporosis. Potential limi- 
tations of GWASs and their meta-analyses are evaluated, with an emphasis on under- 
standing the reasons for inconsistent results between different studies and clarification 
of misinterpretation of GWAS meta-analysis results. Implications and challenges of 
GWAS are also discussed, including the need for multi- and inter-disciplinary studies. 

Keywords: Genome-wide association study, Osteoporosis 



INTRODUCTION 

In the past 4 years, genome-wide association studies (GWASs), assaying hun- 
dreds of thousands of single nucleotide polymorphisms (SNPs) in thousands of 
individuals, have identified a number of promising genetic variants that are asso- 
ciated with osteoporosis and related traits. The first published whole genome se- 
quencing study for osteoporosis and fractures identified a novel rare nonsense 
mutation. This article reviews the current status of GWASs of osteoporosis and 
their meta-analyses with an emphasis on prominent results, approaches, and 
problems with these studies. The major findings of the first whole genome se- 
quencing study for osteoporosis and fractures are also briefly discussed. We focus 
primarily on bone mineral density (BMD), the most important risk trait for osteo- 
porosis, and on osteoporotic fracture (OF), the most severe clinical outcome of os- 
teoporosis. We address how to interpret the discordance of research findings be- 
tween individual GWASs and meta-analyses, and between different meta-analy- 
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ses. The values and limitations of GWAS and meta-analysis 
are evaluated based on empiric and theoretical analyses. 
Finally, future directions using multi- and inter-disciplinary 
study strategies for genetic research of osteoporosis are 
discussed, with potentially significant implications for the 
general human genetics field. 

GWASs AND META-ANALYSIS ON 
OSTEOPOROSIS AND OF 

Osteoporosis is the most common metabolic skeletal dis- 
order in humans. It predisposes people to fragility fractures 
and confers substantial morbidity and mortality, affecting 
over 200 million people worldwide.[1,2] Osteoporosis is 
mainly characterized by low BMD, a highly heritable trait 
with heritability ranging from 0.5 to 0.8.[3-6] OF, as an end- 
point clinical outcome of osteoporosis, also has moderate 
heritability, of approximately 0.5-0.7.[7,8] To date, GWASs 
and their meta-analyses have identified over 60 genes/loci 
associated with variations in BMD and more than 20 genes/ 
loci associated with risk of OF. In addition, a most recently 
published whole-genome sequencing study identified a 
rare nonsense mutation novel within a novel gene LGR4 
that was strongly associated with low BMD and OF.[9] 

The majority of published GWASs have focused on BMD 
using SNP data. A recent review by Richards et al.[10] sum- 
marized the major findings from SNP-based GWASs, which 
will not be repeated here. Instead, we highlight prominent 
genes or loci identified in SNP based GWASs with a focus 
on consistency and inconsistency of results. We address is- 
sues related to interpretation of meta-analysis results and 
replication of study findings among GWASs and meta- 
analyses. A few GWASs have focused on genome wide 
analyses of copy number variants (CNV) and biological 
pathways. The major findings from these latter studies are 
briefly introduced and discussed. 

1. GWASs based on SNPs 

To date, a total of 1 9 GWASs have been published for os- 
teoporosis. Of these, 14 are individual GWASs, and 5 are 
GWAS meta-analyses. The significant genes and loci identi- 
fied from these studies, along with information regarding 
study design (e.g., individual GWASs or meta-analysis) and 
phenotypes (e.g., hip and/or spine BMD or OF) are summa- 
rized in tables 1 and 2. 



1) GWAS meta-analyses 

The Genetic Factors of Osteoporosis (GEFOS) consortium 
published two large-scale GWAS meta-analyses.[11,12] 
Theirfirst meta-analysis (GEFOS-1), which included 19,195 
subjects of European descent, identified 20 BMD loci.[11] 
Their second GWAS meta-analysis (GEFOS-2), the largest 
one to date in the bone field, included 32,961 individuals 
in the discovery phase and was replicated in 50,933 inde- 
pendent subjects.D 2] The study subjects included Europe- 
ans and East Asians. GEFOS-2 identified a total of 56 BMD 
loci at the genome-wide significance level. Of these 56 
loci, 32 were novel and the remaining 24 were genes/loci 
that were previously known to affect bone mass regulation 
and metabolism (e.g., receptor activator of nuclear factor 
(NF)-kappaB (kB) [RANK], RANK ligand [RANKL], and lipo- 
protein receptor-related protein 5 [LRP5]). Notably, multi- 
ple loci achieved highly significant associations in GE- 
FOS-2, with the magnitude of P values reaching10- 60 for 2 
loci, 1 0- 50 for 1 locus, 1 0- 40 for 2 loci, 1 0- 30 for 6 loci, 1 0- 20 for 
1 1 loci, and 10- 10 for 27 loci. GEFOS-2 also revealed 14 loci 
associated with risk of fractures (Table 2). However, frac- 
tures used in the analyses were quite heterogeneous, and 
included hip, spine, and wrist, as well as other types of 
fractures. Due to the well-known genetic and non-genetic 
etiological heterogeneity underlying different types of 
fractures,[2,8,13-16] these study findings should be inter- 
preted with caution, prior to independent validation in 
other samples with homogeneous fracture types. 

Koller et al.[1 7] carried out a meta-analysis of GWASs re- 
stricted to premenopausal white women from four cohorts 
(n =4,061, aged 20 to 45 years), a subset of GEFOS-2, and 
identified two loci (wingless-type MMTV integration site 
family, member 16 [WNT16] and estrogen receptor 1 
[ESRl]IC6orf97] influencing peak bone mass at the lumbar 
spine and femoral neck. Only 4 of the 56 loci detected in 
the joint female GEFOS analysis[1 2] were observed to have 
P values below 5x1 0 s in this study. 

Zhang et al.[18] conducted a three-stage GWAS meta- 
analysis in 27,061 study subjects. Stage 1 meta-analyzed 
seven GWA samples and 1 1,140 subjects for BMDs at the 
lumbar spine, hip and femoral neck, followed by a Stage 2 
in silico replication of 33 SNPs in 9,258 subjects, and by a 
Stage 3 de novo validation of three SNPs in 6,663 subjects. 
Combining evidence from all the stages, two novel loci 
were identified at the genome-wide significance level: 
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Table 1. Putative bone mineral density genes identified in genome-wide association studies and meta-analyses 



Gene 


Full name 


Most significant 
Pvalue 


Initial discovery 
study design 


Reference 


Initial discovery 
population 


Within-study 
replication 


Between-study 
replication 


ABCF2 


ATP-binding cassette, 
sub-family F 


2.56X101-5) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




AD- 

AMTS18 


ADAM metallopeptidase 
with thrombospondin type 
1 motif, 18 


2.56X1 0(-5) 


GWAS 


Xiong etal.[30] 


US Caucasians 


Chinese 


Koller et al.[31] 
in US Caucasians 
and African 
Americans 


ALDH7A1 


Aldehyde dehydrogenase 
7 family, member A1 


2.08x10(-9) 


GWAS 


Guo etal. [37] 


Chinese 


Caucasians and 
Chinese 




ANAPC1 


Anaphase promoting 
complex subunit 1 


1.5x10(-9) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




ARHGAP1 


Rho GTPase activating 
protein 1 


4.0X10I-9) 


Meta-analysis 


Rivadeneira et al. 
[11] 


Caucasians of 
Northern Euro- 
pean 




Estrada etal. [12] 
in Caucasians and 
Asians 


AXIN1 


Axin 1 


2.2x10(-8) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




C7orf58 


Chromosome 7 open 
reading frame 58 


1.1x10(-9) 


Meta-analysis 


Estrada etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




C120RF23 


Chromosome 12 open 
reading frame 23 


9.6X10(-10) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




C180RF19 


Chromosome 18 open 
reading frame 19 


4.9x10(-8) 


Meta-analysis 


Estrada etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




CDKAL1 


CDK5 regulatory subunit 
associated protein 1 -like 1 


1.8x10(-8) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




CLCN7 


Chloride channel, 
voltage-sensitive 7 


1.5X10(-16) 


GWAS 


Duncan et al.[26] 


Postmenopausal 
Caucasian women 


Caucasians 


Estrada etal. [12] 
in Caucasians and 
Asians 


CLDN14 


Claudin 14 


4.15x10(-9) 


Meta-analysis 


Zhang etal.[18] 


Caucasians and 
Asians 


Caucasians and 
Asians 


Estrada etal. [12] 
in Caucasians and 
Asians 


CPN1 


Carboxypeptidase N, 
polypeptide 1 


9X10(-10) 


Meta-analysis 


Estrada etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




CRHR1 


Corticotropin releasing 
hormone receptor 1 


1.4x10(-8) 


Meta-analysis 


Rivadeneira et al. 
[11] 


Caucasians of 
Northern Euro- 
pean 






CTNNB1 


Catenin (cadherin-associ- 
ated protein), beta 1, 


8.1 X10(-10) 


Meta-analysis 


Rivadeneira et al. 
[11] 


Caucasians of 
Northern Euro- 
pean 




Estrada etal. [12] 
in Caucasians and 
Asians 


CYLD 


Cylindromatosis 
(turban tumor syndrome) 


6.2x10(-8) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




DCDC5 


Doublecortin domain 
containing 5 


2.3X101-9) 


Meta-analysis 


Rivadeneira et al. 
[11] 


Caucasians of 
Northern Euro- 
pean 




Estrada etal. [12] 
in Caucasians and 
Asians 


DHH 


Desert hedgehog 


1.2X 10(-15) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




DNM3 


Dynamin 3 


8.5X10(-15) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




D0K6 


Docking protein 6 


8.87x10(-7) 


GWAS 


Hsu etal.[27] 


US Caucasians 


Caucasians of 

Northern 

European 




ERC1 


ELKS/RAB6 interacting/ 
CAST family member 1 


5.6X10(-12) 


Meta-analysis 


Estrada etal. [12] 


Caucasians and 
Asians 


Caucasians and 
Asians 
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Table 1. (Continued from the previous page) Putative bone mineral density genes identified in genome-wide association studies and meta-analyses 



Gene 


Full name 


Most significant Initial discovery 
Pvalue study design 


Reference 


Initial HKrnvpp./ 

ii ii li d i uioijUvci y 

population 


Within-sti iHv 

vviuiiii oLuuy 

replication 


Rptwppn-sti iHv 

UGLVVGGII oLUUy 

replication 


ESR1 


Estrogen receptor 1 


6.1x10(-11) 


GWAS 


Styrkarsdottir 

at al r?/ll 
gL dl.[Z4J 


Icelandic 

l^ai m a c i a n c 
LdULdbldTlb 


Caucasians of 

Wcblclll EUlUpc. 


Rivadeneira et al.[1 1] 

and RinharHc at al I9Q1 
dllU MlLlldrUb eL dl.[ZJJ 

in Caucasians. Koller 
et al.[31] in premeno- 
pausal women. 


FAM9B 


Family with sequence 
similarity 9, member B 


1.2x10(-8) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




FLJ42280 


Putative uncharacterized 

pi ULcll I 


3.4X10(-12) 


Meta-analysis 


Rivadeneira 

at al M 11 
GL dl .[ 1 1 J 


Caucasians of 

ft nrthorn Piirn 
IM UI U lei II CUI U 

pean 




Estrada etal.[12] in 

UdULdbldllb dllU Hbldllb 


nvtiMZ/ 
GREM2 


rurmin z/ gremlin z 


1.9x10(-9) 


GWAS 


rdLcmosier 
etal.[19] 


UdULdbldllb 01 

Northern Euro- 
pean 


l^aimaoianc 
UdULdbldllb 01 

Northern Euro- 
pean 




FOXL1 


Forkhead box L1 


1.7x10(-8) 


Meta-analysis 


Rivadeneira 

at al [1 11 
GL dl .[ 1 1 J 


Caucasians of 

K nrtnorn E-iiiyi 
IM UI U lcl II CUI U 

pean 




Styrkarsdottir etal.[70] 

in Pact Acianc PctraHa 
III CdbL Hbldllb. CbLldUd 

et al.[12] in Caucasians 
and Asians 


FUBP3 


Far upstream element (FUSE) 
binding protein 3 


3.1 X10(-8| 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




GALNT3 


UDP-N-acetyl-alpha-D- 
galactosamine:polypeptide 
N-acetylgalactosaminyl- 
transferase 3 


2.3X101-5) 


GWAS 


Duncan 
et al.[26] 


Postmenopausal 
Caucasian women 


Caucasians 


Estrada et al.[1 2] in 
Caucasians and Asians 


GPATCH1 


G patch domain containing 1 


9.2x10(-9) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




GPR177/ 


WLS 


Wntless homolog 
(Drosophila) 


3.3x10(-13) 


Meta-analysis 


Rivadeneira 
etal.[11] 


Caucasians of 
Northern Euro- 
pean 




Styrkarsdottir et al.[70] 
in East Asians. Hsu et 
al.[27] in US Cauca- 
sians. Duncan et al.[26] 
in Postmenopausal 
Caucasian women. 
Estrada etal.[12] in 
Caucasians and Asians. 


HDAC5 


Histone deacetylase 5 


1.7x10(-8) 


Meta-analysis 


Rivadeneira 
etal.[11] 


Caucasians of 
Northern Euro- 
pean 






IBSP 


Integrin-binding sialoprotein 


7.6x10(-7) 


GWAS 


Duncan 
etal.[26] 


Postmenopausal 
Caucasian wom- 
en 


Caucasians 




IDUA 




5.2X10(-15) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




IL21R 


Interleukin 21 receptor 


2.31 X10(-B) 


GWAS 


Guo etal.[37] 


Chinese 


Caucasians and 
Chinese 




INSIG2 


Alpha L-iduronidase 


1.2x10(-10) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




JAG1 


Jagged 1 


5.27x10(-8) 


GWAS 


Kung et al.[71] Southern Chinese 
women 




Estrada et al.[1 2] in 
Caucasians and Asians 


KIAA2018 




4.1 X10(-10) 




Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




KCNMA1 


Potassium large conduc- 
tance calcium-activated 
channel, subfamily M, alpha 
member 1 


1.5x10(-12) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 
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Table 1. (Continued from the previous page) Putative bone mineral density genes identified in genome-wide association studies and meta-analyses 



Gene 


Full name 


Most significant Initial discovery 
Pvalue study design 


Reference 


Initial discovery 
population 


Within-study 
replication 


Between-study 
replication 


LACTB2 


Lactamase, beta 2 


4.7XlO(-8) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




LEKR1 


Leucine, glutamate and ly- 
sine rich 1 


4.5x 10(-12) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




LIN7C 


lin 7 homologue C 


4.9x10(-8) 


Meta-analysis 


Estrada 

pt al [1 91 
cL dl.[ 1 L\ 


Caucasians and 

Mbldl lb 


Caucasians and 
Asians 




LRP4 


Low density lipoprotein 
receptor-related protein 4 


4.0x10(-9) 


Meta-analysis 


Rivadeneira 
etal.[11] 


Caucasians of 
Northern Euro- 
pean 




Estrada etal. [12] in 
Caucasians and Asians 


LRP5 


Low density lipoprotein 
receptor-related protein 5 


6.3X10(-12) 


GWAS 


Richards 
etal.[23] 


Caucasians 


Caucasians of Rivadeneira et al.[1 1] in 
western Europe Caucasians of Northern 
European. Styrkarsdot- 
tir etal. [70] in East- 
Asians 


LRRC4C 


Leucine rich repeat contain- 

inn AC 

iny ^ii 


8.89x10(-7) 


GWAS 


Hsu etal.[27] 


US Caucasians 


Caucasians of 
Northern Euro- 
pean 




MARK3 


MAP/microtubule affinity- 
regulating kinase 3 


1.8x10(-9) 


GWAS 


Styrkarsdottir 
etal.[72] 


Icelanders 


Caucasians of 
European de- 
scent 


Rivadeneira etal. [11] in 
Caucasians of Northern 
European. Estrada et 
al.[12] in Caucasians 
and Asians 


MBL2 


Mannose-binding lectin 

nrntoin fi *7 cnlnHIo 
lUIULclll LiJ L, oUIUUIc 


1.BX10(-12) 


Meta-analysis 


Estrada 

pt al M 91 
cL dl.[ 1 L\ 


Caucasians and 

Acionc 


Caucasians and 
Asians 




MEF2C 


Myocyte enhancer factor 2C 


1.3x10(-13) 


GWAS 


Duncan 
et al.[26] 


Postmenopausal 
Caucasian wom- 
en 


Caucasians 


Rivadeneira etal. [11] in 
Caucasians of Northern 
European. Estrada et 
al.[12] in Caucasians 
and Asians 


IVILrL 


ividxrix exirdLciiuidr pnus- 
phoglycoprotein 


4.0 X 10(-9) 


Meta-analysis 


nlvdUcllcird 

etal.[11] 


LdULdSldHS 01 

Northern Euro- 
pean. 




Styrkarsdottir et al.[70] 
in East-Asians 


MPP7 


Membrane protein, palmi- 
toylated 7 (MAGUK p55 
subfamily member 7) 


2.9X 10(-9) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




0SBPL1A 


Oxysterol binding protein- 
like 1 A 


Ml X 10(-/) 


bWAb 


Hsu etal.[27] 


US Caucasians 


Caucasians of 
Northern Euro- 
pean 




PTHLH 


Parathyroid hormone-like 
hormone 


1.9X10(-12) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




RAP1A 


RAS-related protein RAP1A 


2.80 X 10(-7) 


r*i MAO 

GWAS 


Hsu etal.[27] 


US Caucasians 


Caucasians of 
Northern Euro- 
pean 




PKDCC 


Protein kinase domain con- 
taining, cytoplasmic homo- 
logue 


1.3x10(-9) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




NTAN1 


Ribosomal protein S6 kinase, 
90kda, polypeptide 5 


1.7X101-11) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




RSP03 


R-spondin 3 


2.2X101-9) 


GWAS 


Duncan 
et al.[26] 


Postmenopausal 
Caucasian wom- 
en 


Caucasians 


Estrada etal.[12] in 
Caucasians and Asians 


RUNX2 


Runt-related transcription 
factor 2 


5.6x10(-11) 


Meta-analysis 


Caucasians 
and Asians 


Caucasians and 
Asians 


Caucasians and 
Asians 





(Continued to the next page) 
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Table 1. (Continued from the previous page) Putative bone mineral density genes identified in genome-wide association studies and meta-analyses 



Gene 


Full name 


Most significant Initial discovery 
Pvalue study design 


Reference 


Initial discovery 
population 


Within-study 
replication 


Between-study 
replication 


SALL1/ 
CYLD 


Sal-like 1 (Drosophila) 


5.0X10(-12) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




SLC25A13 


Solute carrier family 25 


a.l X 1U(-48) 


Meta-analysis 


Estrada 
etal.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




OIVIuo 


oi i iy-D iiomuioy, nuiiociioc 
mediated mrna decay factor 
(C. Elegans) 


1.7X101-8) 


Meta-analysis 


CbLldUd 

etal.[12] 


UdULdbldilb dllU 

Asians 


UdULdbldilb dill] 

Asians 




OIVIUu 1 


QPARP rolatoH mnHiilar 

orMnb reidieu lilUUUIdr 

calcium binding 1 


3.98X10(-13) 


Meta-analysis 


Zhang 
etal.[18] 


UdULdbldilb dllU 

Asians 


UdULdbldilb dllU 

Asians 


FetraHa at al M 91 in 
LbudUd el dl.[ 1 ZJ 111 

Caucasians and Asians 


SOST 


Sclerostin 


2.1X101-8) 


GWAS 


Styrkarsdottir 
etal.[72] 


Icelanders 


Caucasians of 
European de- 
scent 


Styrkarsdottir etal.[70] 
in East-Asians. Estrada 

at qI r 1 *7 1 in Pqi ip^cionc 
cl dl-L 1 L\ III UdULdbldilb 

and Asians 


OUAH 


QRV loov nQtormininn romnn 

on I \bcA ucici 1 1 111 111 iy i uy iui i 

Y)-box4 


1.8x10(-8) 


GWAS 


UUI ILdl 1 

etal.[26] 


Pnotmannnai loa 1 
rUbLlllGllUUdUbdl 

Caucasian wom- 
en 


UdULdbldl lb 


FetraHa Pt al M 71 in 
CbU dUd c L dl .[ 1 L\ III 

Caucasians and Asians. 


S0X6 


SRY (sex determining region 

Yl hnv R 


6.4x10(-10) 


Meta-analysis 


Rivadeneira 

Pt al Mil 
c I d 1 .[ 1 1 J 


Caucasians of 

N nrtnQrn f-iirn 
INIUI 11 lcl 1 1 CUI U 

pean 




Estrada etal.[12] in 

UdULdbldilb dllU Hbldllb. 

Liu et al.[48] in US Cau- 
casians 


OUAj 


on i-uux Luiudiniiig gene y 


3.7x10(-8) 


Meta-analysis 


CbLldUd 

etal.[12] 


UdULdbldilb dllU 

Asians. 


UdULdbldilb dllU 

Asians. 




SP7 


Sp7 transcription factor 7 


9.9x10(-11) 


GWAS 


Styrkarsdottir 

Pt al T711 
cL al. |7 1 J 


Icelanders 


Caucasians of 

Fi irnnQon Ho 
CUIUptrdll Uc 

scent 


Timpson et al.[32] in 

[ hilHran RiwaHanaira 
Ul II IUI cl 1. niVdUcllclld 

et al.[1 1] in Caucasians 
of Northern European 


orr 1 


Secreted phosphoprotein 1 


6.0x10( -8) 


GWAS 


111 in/^on 

UUllLdll 

etal.[26] 


r (JSLlllcll(J|JdUbdl 

Caucasian wom- 
en 


UdULdbldilb 




orr z 


Secreted phosphoprotein 2, 4.64 x 1 0(-7) 
24 kda 


GWAS 


Uqii pt q| [971 

nbU cL a\.[L/ J 


1 IQ roi ipccionc 
UO LidULdbldl lb 


Panpacianc nT 
LidULdbldl lb Ul 

Northern Euro- 
pean 




or i d/v / 


Spectrin, beta, non-erythro- 
cytic 1 


1.6x10(-8) 


Meta-analysis 


R i wQrlonQi ro 
Ml VdLlcl Icll d 

etal.[11] 


Paiipaciano nT 
UdULdbldl lb U 

Northern Euro- 
pean 




FetraHa at al M?l in 
Lb LI dUd G L dl .[ I ZJ III 

Caucasians and Asians. 


STARD3NL 


STARD3 N-terminal like 


1.1X101-9) 


Meta-analysis 


Rivadeneira 
etal.[11] 


Caucasians of 
Northern Euro- 
pean 




Estrada etal.[12] in 
Caucasians and Asians 


TBC1D8 


TBC1 domain family, mem- 
ber 8 (with GRAM domain) 


1.48X101-7) 


GWAS 


Hsu etal.[27] 


US Caucasians 


Caucasians of 
Northern Euro- 
pean 




TGFBR3 


Transforming growth factor, 2.13 x 1 0(-8) 
beta receptor III 


GWAS 


Xiong 
etal.[30] 


US Caucasians 


Chinese 




TNFSF1 1/ 
RANKL 


Tumor necrosis factor 
(ligand) superfamily, mem- 
ber 11 


2.0X10(-14) 


GWAS 


Richards 
etal.[23] 


Caucasians 


Caucasians of Duncan et al. [26] 
western Europe in Postmenopausal 
Caucasian women. 
Styrkarsdottir et al.[24] 
in Caucasians. Estrada 
et al.[12] in Caucasians 
and Asians. 



(Continued to the next page) 
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Table 1. (Continued from the previous page) Putative bone mineral density genes identified in genome-wide association studies and meta-analyses 



Gene 


Full name 


Most significant Initial discovery 
Pvalue study design 


Reference 


Initial discovery 
population 


Within-study 
replication 


Between-study 
replication 


TNFRS- 

F11A/ 

RANK 


Tumor necrosis factor recep- 
tor superfamily, member 
11a, NFKB activator 


9.4x10(-9) 


GWAS 


Styrkarsdottir 
eta I. [24] 


Icelandic 
Caucasians 


Caucasians of 
European de- 
scent 


Richards etal.[23] in 
Caucasians of western 
Europe. Rivadeneira et 
al.pl] in Caucasians 
of Northern European. 
Estrada etal.[12] in 
Caucasians and Asians. 


TNFRS- 

F11B/ 

ORG 


Tumor necrosis factor recep- 
tor superfamily, member lib 


3.5x10(-16) 


GWAS 


Styrkarsdottir 
etal.[24] 


Icelandic 
Caucasians 


Caucasians of 
European de- 
scent 


Richards etal.[23] in 
Caucasians of western 
Europe. Rivadeneira et 
al.P1] in Caucasians 
of Northern European. 
Estrada etal.[12] in 
Caucasians and Asians. 


WNT16 


Wingless-type MMTV inte- 
gration site family, member 
16 


1.4x10(-16) 


Meta-analysis 


Estrada et 
al.[12] 


Caucasians and 
Asians 


Caucasians and Koller et al.[17] in pre- 
Asians menopausal women 


WNT4 


Wingless-type MMTV inte- 
gration site family, member 
4 


9.6X10(-11) 


Meta-analysis 


Rivadeneira 
etal.[11] 


Caucasians of 
Northern Euro- 
pean 




Duncan et al.[26] 
in postmenopausal 
Caucasian women. 
Styrkarsdottir et al.[70] 
in East-Asians. Estrada 
et al.[12] in Caucasians 
and Asians. 


XKR9 


XK, Kell blood group com- 
plex subunit-related family, 
member 9 


4.7X101-8) 


Meta-analysis 


Estrada et 
al.[12] 


Caucasians and 
Asians 


Caucasians and 
Asians 




ZBTB40 


Zinc finger and BTB domain 
containing 40 


3.2x10(-10] 


Meta-analysis 


Rivadeneira 
etal.[11] 


Caucasians of 
Northern Euro- 
pean 




Duncan et al.[26] 
in postmenopausal 
Caucasian women. 
Styrkarsdottir et al.[70] 
in East-Asians. Estrada 
et al.[12] in Caucasians 
and Asians. 



Genes/loci in bold are those with evidence of cross study replication. 

Within-study replication means the study includes a follow-up replication component. 

Between-study replication means the gene/loci showed significance in different studies (not necessarily totally independent because some large meta- 
analyses such as GEFOS-2 (Estrada et al. (1 2)) include the samples from other individual GWASs. 
LS, lumbar spine; FN, femoral neck; OF, osteoporotic fracture; HF, hip fracture. 



14q24.2 (rs227425, P=3.98x10 13 , SPARC related modular 
calcium binding 1 [SMOG]) and 21 q22.1 3 (rs170183, 
P=4.15x10" 9 , claudin 14 [CLDN14]). These two SNPs were 
also significant in GEFOS-2.[12] This study independently 
confirmed 13 previously reported loci. [18] Further gene 
expression analysis in osteogenic cells implied potential 
functional association of the two novel candidate genes 
(SMOG and CLDN14) in bone metabolism. 

Most studies have focused on areal BMD (aBMD) ob- 
tained from a 2-dimensional projectional scan with dual 
energy X-ray absorptiometry (DXA). Although aBMD is the 
gold standard for diagnosing osteoporosis, it fails to pro- 
vide a detailed information necessary to discern traits such 



as trabecular volumetric BMD (vBMD), cortical vBMD and 
bone microstructural parameters. Quantitative computed 
tomography (QCT) analysis has the advantage to reveal 
unique information about these bone traits. Paternoster et 
al.[19] published the first GWAS to identify genetic loci as- 
sociated with cortical and trabecular bone microstructural 
parameters in European Caucasians. Their cortical vBMD 
GWAS meta-analysis (n = 5,878) followed by replication (n 
= 1,052) identified genetic variants in four separate loci 
(RANKL, rs1021 188, P = 3.6x10 14 ; LOC285735, rs271 170, 
P=2.7x10 12 ; osteoprotegerin [OPG], rs7839059, P=1.2 x 
10" 10 ; and ESR1/C6orf97, rs6909279, P=1.1x10" 9 ).The tra- 
becular vBMD GWA meta-analysis (n =2,500) followed by 
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Table 2. Putative osteoporotic fracture genes identified in genome-wide association studies and meta-analyses 



Gene 


Full name 


Phenotype 


oiuuy uesign 
(GWAS, meta-analysis) 


Reference 


[VIOSl SigillllLdlU 

Pvalue 


ADAMTS18 


ADAM metallopeptidase with thrombospondin 
type 1 motif, 18 


Hip-OF 


GWAS 


Xiong et al 


[30] 


2.9 X 101-2) 


ALDH7A1 


Aldehyde dehydrogenase 7 family, member A1 


Hip-OF 


GWAS 


Guo et al 


[37] 


8.53x10(-9) 


C17orf53 


Chromosome 17 open reading frame 53 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


4.1X10I-5) 


CTNNB1 


Catenin (cadherin-associated protein), beta 1 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


2.9X1 0(-7) 


DCDC5 


Doublecortin domain containing 5 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


3.3X10I-5) 


FAM2WA 


Family with sequence similarity 210, member a 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


8.8X10(-13) 


FUBP3 


Far upstream element (FUSE) binding protein 3 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


3.5X101-5) 


LRP5 


Low density lipoprotein receptor-related protein 5 Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


1 .4XHK-8) 


MBL2/DKK1 


Lectin, mannose-binding, soluble, 2 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


9.0x10(-9) 


MECOM 


MDS1 and EVI1 complex locus 


OF 


GWAS 


Hwang et al 


[38] 


3.59x10(-8) 


MEPE/SPP1 


Matrix, extracellular, phosphoglycoprotein 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


1.7x10(-8) 


RPSBKA5 


Ribosomal protein S6 kinase, polypeptide 5 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


7.2X1 0(-5) 


SLC25A13 


Solute carrier family 25 (citrin), member 1 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


5.9x10(-11) 


CDCT 


Sclerostin 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


c n v 1 n/ c\ 
b.y x i u(-o] 


QDTDM1 

or I UN I 


spectrin, beta, nonerythrocytic, 1 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


Z.bX( IU(-o] 


QTAPFV3MI 


STARD3 N-terminal like 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


1 .1 X 1 U(-Dj 


/ brbnj 


Transforming growth factor, beta receptor III 


Hip-OF 


bVVAo 


Xiong et al 


[30] 


L. I o X I U(-oJ 


TNFRSF1 1A 
(RANK) 


Tumor necrosis factor receptor superfamily, 
member 1 1 a, NFKB activator 


NV/V-Fracture 


Meta-analysis 


Richards et al 


[23] 


0.02 


UGT2B17 


UDP glucuronosyltransferase 2 family, 
polypeptide B17 


Hip OF 


GWAS 


Yang et al 


[42] 


5.0x10(-4) 


WNT4 


Wingless-related MMTV integration site 4 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


1.4x10(-7) 


WNTW 


Wingless-related MMTV integration site 16 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


2.7x10(-7) 


ZBTB40 


Zinc finger and BTB domain containing 40 


Any type of fracture 


Meta-analysis 


Estrada et al 


[12] 


3.6x10(-6) 



Any type of fracture: consisting of low-trauma fractures at any skeletal site (except fingers, toes and skull). 
NV Fracture, Nonvertebral fracture; OF, Osteoporotic fracture; GWAS, genome-wide association study. 



replication (n = 1,022) identified one locus reaching ge- 
nome-wide significance (formin 2 [FMN2]/ gremlin 2, DAN 
family BMP antagonist [GREM2], rs9287237, P=1.9x10" 9 ). 
In addition, rsl 021 1 88 was associated with cortical porosi- 
ty while rs9287237 was associated with trabecular bone 
fraction. The genetic variant in the FMN2/GREM2 locus was 
also associated with fracture risk in the MrOS Sweden co- 
hort and GREM2 expression in human osteoblasts. Two of 
these (FMN2/GREM2 and LOC285735) are novel bone re- 
lated loci, while the other three have previously been re- 
ported to be associated with aBMD. This study provided 
evidence that the genetic determinants of cortical and tra- 
becular vBMDs differ. However, QCT has its limitations, in- 
cluding being not applicable to World Health Organization 
(WHO) definition of osteoporosis that is based on DXA mea- 
surement, being more expensive with a higher dosage of 
exposure to radiation and may not predict fractures better 



than DXA measurement.[20,21] Nevertheless, its advan- 
tages over DXA make QCT a complementary (but not nec- 
essarily replacement) approach to bone health assess- 
ment.^] 

Conducting these large-scale meta-analyses requires an 
extensive, concerted effort in collaboration and coordina- 
tion among various research centers and groups, as well as 
centralization and standardization in data analyses. Thus, 
publication of these studies represents one of the most 
impressive achievements in the osteoporosis genetics 
field. Based primarily on the large samples involved in me- 
ta-analyses such as these, some researchers have begun to 
consider large meta-analyses as a gold standard for evalu- 
ating other individual GWASs, with an explicit assumption 
that validity of the findings from other studies of smaller 
samples needs to be confirmed in large meta-analyses 
such asGEFOS-2.[12] 
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Comparing the results of the two published GEFOS me- 
ta-analyses, it is interesting to note that several loci identi- 
fied in GEFOS-1 (e.g., corticotropin releasing hormone re- 
ceptor 1 [CRHR1] and histone deacetylase 5 [HDAC5]) were 
not significant in GEFOS-2, despite the fact that the num- 
ber of samples of GEFOS-2 almost tripled the number in 
GEFOS-1, and the majority of samples used in GEFOS-1 
were included in GEFOS-2. 

Further comparing the results of these meta-analyses 
with individual GWASs, it is intriguing that some significant 
loci identified in meta-analyses were found to be signifi- 
cant at the genome-wide level in individual GWASs, while 
others were not. Likewise, some significant loci identified 
in individual GWASs were replicated in meta-analyses, 
while many others were not. We will highlight some of 
these consistencies and inconsistencies in the following 
paragraphs and discuss how to interpret these findings. 

2) Individual GWASs 

(1) Single ethnicity 

The first two individual GWASs were performed in hu- 
man subjects of European ancestry [23,24] using the data 
from TwinsUK/Rotterdam and deCODE Genetic studies. 
These two studies identified a total of five loci at the ge- 
nome-wide significance level (P< 5x1 0 s ). The associated 
genes, including LRP5, OPG, RANK, RANKL and ESR1, are 
well known for their importance to bone and mineral me- 
tabolism, and have been well established in previous can- 
didate gene studies. 

The efficiency of extreme-truncate selection design for 
quantitative trait association studies has been well estab- 
lished.^] Duncan et al. adopted this design by focusing 
on postmenopausal Caucasian women with either ex- 
tremely high or low hip BMD in the discovery stage.[26] Al- 
though the sample size of the discovery phase was moder- 
ate (1,955 subjects in discovery phase), the power of the 
study was still high due to the extremely discordant sample 
design. In addition to replicating 21 of 26 known BMD-as- 
sociated genes, they identified six new associations in/ 
around the genes chloride channel, voltage-sensitive 7 
(CLCN7), UDP-N-acetyl-alpha-D-galactosamine:polypeptide 
N-acetylgalactosaminyltransferase 3 (GALNT3), integrin- 
binding sialoprotein {IBSP), latent transforming growth fac- 
tor beta binding protein 3 (LTBP3), R-spondin 3 (RSP03), 
and SRY (sex determining region Y)-box 4 (50X4). Some of 



these genes (e.g., LTBP3) were not significant at the ge- 
nome-wide significance level in GEFOS-2, although the 
samples used in this study were included in GEFOS-2. 

Organic integration of information from different levels 
and disciplines may help identify genetic variants of mod- 
est or moderate effects. By integrating GWAS in Caucasians 
and expression signature profiling of various human tis- 
sues, Hsu et al.[27] identified three novel genes, member 
of RAS oncogene family (RAP1A), TBC1 domain family, 
member 8 (with GRAM domain) {TBC1D8), and oxysterol 
binding protein-like M\(0SBPL1A), for BMD, and replicated 
a well-known gene, OPG. They also prioritized 16 sugges- 
tive genome-wide significant candidate genes, among 
which two candidate genes, G protein-coupled receptor 
177 (GPR177) and SOX6, were replicated in GEFOS-2. The 
samples used in this study (mainly from the Framingham 
Osteoporosis Study) were included in GEFOS-1 and GE- 
FOS-2. Therefore, it does not represent an independent 
replication. The three identified novel genes (i.e., RAPT A, 
TBC1D8, and OSBPL1A) were not significant at the ge- 
nome-wide level in either GEFOS-1 [11] or GEFOS-2.[1 2] 
Gene expression profiling analyses in this study were con- 
ducted in lymphocytes and the liver whose expression 
profiles may not be directly related to bone. Moreover, it is 
the differential gene expression with regard to osteoporo- 
sis related phenotypes, not simply expression of genes, 
which is more clinically relevant to risk of osteoporosis and 
OF.[28,29] 

(2) Multi-ethnicities 

To investigate potential population generality/specificity 
of genes/loci, several GWASs were conducted in popula- 
tions of multiple ethnicities. Xiong et al.[30] used a US 
Caucasian cohort as the discovery sample followed by rep- 
lication analyses in independent populations including 
those of African (from the Tobago study) and Asian ances- 
try. Two novel genes, ADAM metallopeptidase with throm- 
bospondin type 1 motif, 1 8 (ADAMTS1 8) and transforming 
growth factor, beta receptor III (TGFBR3), were identified at 
the genome-wide significance level. These two genes were 
found to be significant in other individual GWASs conduct- 
ed in completely independent samples by independent 
investigators in premenopausal US Caucasian women[31] 
and in postmenopausal women with extremely discordant 
BMD values (28). However, neither of these genes was 
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found to be significant in the two GEFOS meta-analyses. 
[11,12] 

Focusing on subjects of a specific gender and/or age 
range may reduce the potential confounding effects of 
heterogeneity, and empirically increase the power of a 
study. Koller et al.[31] conducted a GWAS in 1,524 pre- 
menopausal US Caucasian women aged 20-45 years and 
669 African American premenopausal women aged 20-44 
years. A novel gene, catsper channel auxiliary subunit beta 
{CATSPERB), was identified as being significant for peak 
BMD. Despite the fact that samples from this study were 
included in GEFOS-2,[12] CATSPERB was not found to be 
significant in the meta-analyses. 

(3) Children 

Although the majority of GWASs have been focused ex- 
clusively on adults, Timpson et al.[32] published the first 
GWAS in 1 ,51 8 children from the Avon Longitudinal Study 
of Parents and Children. They identified association with 
BMD at the Osterix (SP7) gene locus, a transcription factor 
responsible for regulating osteoblast differentiation. In an- 
other study, Medina-Gomez et al.[33] identified a novel 
gene, Wntl6 that was associated with total body and skull 
BMD variation in children and adults, which suggested a 
role in peak bone mass accrual which may impact the risk 
of osteoporosis later in life. Interestingly, this locus was 
also associated with cortical bone thickness, wrist BMD, 
bone strength, and risk of forearm fracture in adults,[34] 
peak bone mass in premenopausal women,[17] and BMD 
and fracture in elderly individuals of European descent.D 2] 

(4) OF 

The majority of the published GWASs of osteoporosis 
have been focused on BMD, a major risk factor for OF. How- 
ever, OF is an important clinical end-point of osteoporosis, 
and is the most clinically relevant trait. Earlier studies have 
demonstrated a high genetic determination for OF risk,[1 3] 
and have also shown that BMD and fracture may have dif- 
ferent genetic determination,[7,8] indicating the necessity 
to perform genetic studies on OF perse. Over the past sev- 
eral years, OF, as an independent study phenotype, has re- 
ceived increased attention in the bone genetics field. 
[35,36] Currently, a major problem facing genetic studies 
of OF is the difficulty in recruiting sufficiently large and ho- 
mogeneous samples of the same OF type. 



The first published GWAS for OF used 700 elderly Chi- 
nese Han subjects (350 with hip OF and 350 healthy matched 
controls) in the discovery phase, followed by replication in 
an independent Chinese sample containing 390 cases with 
hip OF and 516 controls.[37] A novel gene, aldehyde dehy- 
drogenase 7 family, member A1 (ALDH7A 7), was identified 
to be significant for hip OF. Further analyses confirmed its 
relevance to hip BMD in both Chinese (sample size, 2,955) 
and Caucasian (sample size, 7,007) populations. Although 
the discovery sample size was relatively small, the study 
was conducted in a homogeneous sample of the same 
ethnicity and the same type of carefully clinically adjudi- 
cated OF, which may empirically increase the power of 
gene identification. Notably, ALDH7A1 was not found to 
be significant in GEFOS-2. 

Hwang et al.[38] published a GWAS for OF in East Asian 
populations. They first conducted a GWAS in the discovery 
cohort which included 288 cases and 1,139 controls, fol- 
lowed by a two-stage de novo replication analyses which 
included 462 cases and 1,745 controls in stage 1 and 369 
cases and 560 controls in stage 2. A new locus associated 
with OF (rs784288 in the MECOM gene) showed genome- 
wide significance (/>= 3.59 X 10 s ; OR 1.39). Further RNA in- 
terference analysis revealed that a MDS1 and EVI1 complex 
locus (MECOM) knockdown suppresses osteoclastogenesis. 

An alternative approach for studying the genetic deter- 
mination of OF is to first identify BMD loci, and then test 
the relevance of those loci to OF. This approach has been 
used for individual GWASs[23,24,30] and for GEFOS-2 (Ta- 
ble 2).[1 2] Using this approach, GEFOS-2 identified 14 BMD 
loci that were associated with risk of OF. One potential 
problem with the GEFOS-2 OF analyses is that it combined 
OFs from a variety of clinically heterogeneous types of OFs, 
thus departing from the assumption of homogeneity 
made by meta-analysis for the phenotype. From perspec- 
tive of bone biology, the proportions of cortical and can- 
cellous bone differ at the different sites in the skeleton 
(e.g., hip, spine and wrist) where OFs frequently occur.[39] 
The differences in behavior of bone at the different sites 
are most likely caused by the different environments of the 
bone cells in cortical or cancellous bone.[39] Moreover, 
each type of fracture has its own unique risk factors.[2] For 
instance, fall is a major risk factor for hip and wrist frac- 
tures; however, it contributes little to the risk of vertebral 
fracture. Epidemiology data also show that the prevalence 
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of fracture at different sites is different in terms of age, 
gender and ethnicity.[2,40] So far, the majority of the ge- 
netic variants identified by GWASs or candidate gene stud- 
ies for BMD at different skeletal sites do not overlap.[5,10] 
Given these evidences, it is likely that different skeletal 
sites may not share much common genetic component. 
Therefore, the mixing of phenotypes in GEFOS-2 meta- 
analysis complicates the interpretation of results, as in- 
creasing sizes of non-homogeneous samples may lead to 
false positive/negative results. Although meta-analysis can 
be very powerful, its application must be executed with 
caution. 

2. GWASs based on CNVs 

SNP-based GWASs, including GEFOS meta-analyses, 
have discovered genes/loci which, collectively, account for 
less than 6% of the risk for developing osteoporosis. There- 
fore, it is important to explore the possibility that some of 
the remaining undiscovered genetic factors that influence 
risk of osteoporosis are due to genomic mechanisms other 
than individual nucleotide changes, such as CNVs. 

CNV is a common type of genomic variability, with varia- 
tions in the size of DNA fragments ranging from 1 kilobase 
to several megabases. CNVs, which may include duplica- 
tions or deletions, can influence gene expression by dis- 
rupting coding sequences, perturbing long-range gene 
regulation, or altering gene dosage, and these effects could 
contribute to phenotypic variations or disease risk.[41 ] 

A CNV-based GWAS for osteoporosis was conducted by 
Yang et al.[42] in 700 elderly Chinese individuals compris- 
ing 350 cases with homogeneous hip OF and 350 matched 
controls. In that study, CNV 4q1 3.2 was strongly associated 
with OF. Further validation experiments identified a dele- 
tion variant of UDP glucuronosyltransferase 2 family, poly- 
peptide B1 7 (UGT2B17) in CNV 4q1 3.2, and the association 
between CNV of UGT2B 77 and OF was replicated in an in- 
dependent Chinese sample containing 399 cases with hip 
OF and 400 controls.[42] Despite the relatively small sam- 
ple sizes, the CNV focused GWAS identified novel genes/ 
loci that were otherwise missed in SNP-based GWASs. 

Oei et al.[43] performed a genome-wide CNV analysis in 
5,178 individuals from Netherlands, including 809 OF cas- 
es, and followed by in silico lookups and de novo replica- 
tion in several independent samples. A rare (population 
prevalence 0.14%) 210 kb deletion located on chromo- 



some 6p25.1 was associated with the risk of fracture. The 
prevalence of this deletion showed geographic diversity, 
being absent in samples from Australia, Canada, Poland, 
Iceland, Denmark, and Sweden, but present in the Nether- 
lands (0.34%), Spain (0.33%), USA (0.23%), England 
(0.15%), Scotland (0.10%), and Ireland (0.06%). Larger and 
geographically restricted studies are needed to confirm 
this regional association. 

3. GWASs based on biological pathways 

Conventional GWASs consider the effects of individual 
genetic markers independently, and focus only on genes 
or markers of top-ranking statistical significance. However, 
genes and/or their products often work together, interact- 
ing in functional groups or pathways to contribute to phe- 
notypic variation or susceptibility to disease. Therefore, 
GWASs which focus on individual genes/loci may not be 
optimally effective for identifying pathophysiological^ 
significant genetic pathways underlying osteoporosis and 
OF. A major advantage of pathway-based GWASs is that 
they link the wealth of information embedded in GWAS 
data to our knowledge of functional biological pathways 
that is readily available in public databases. 

Zhang et al.[44] performed the first pathway-based 
GWAS for osteoporosis in a population of US Caucasians, 
and identified the importance of the regulation-of-au- 
tophagy (ROA) pathway for ultradistal radius BMD. This is 
the first study implicating the ROA pathway's importance 
to osteoporosis, though autophagy had previously been 
recognized to be important in endochondral ossification. 
[45] Another pathway-based GWAS conducted in 700 el- 
derly Chinese Han subjects (350 with hip OF and 350 
healthy matched controls) revealed the importance of the 
tumor necrosis factor (TNF)-related apoptosis-inducing li- 
gand (TRAIL) pathway for hip OF.[46] In the pathway- 
based analyses, the maximum statistic for all SNPs near a 
gene was taken to represent the significance of the gene 
and use a permutation-based approach that shuffles the 
phenotypes to adjust for multiple testing.[47] Notably, 
most individual genes included in these significant path- 
ways did not reach genome-wide significance in SNP- 
based GWASs, underscoring the importance of looking be- 
yond the top significant SNPs/genes and searching for ad- 
ditional risk factors with moderate statistical significance. 
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4. Bivariate GWASs 

Osteoporosis is highly correlated, both genetically and 
environmentally, with other human diseases/conditions 
such as obesity (characterized by high body mass index 
[BMI]) and sarcopenia (characterized by low lean mass and 
strength). GWASs, however, have typically been imple- 
mented in a univariate framework, analyzing different 
phenotypes (e.g., BMD or BMI) independently, and ignor- 
ing the potential genetic correlation between disease 
traits. Thus, this approach is not well suited for detecting 
pleiotropic genes, which may exist for two genetically cor- 
related diseases. One potentially effective strategy for 
identifying pleiotropic genes is to analyze potentially cor- 
related disease phenotypes simultaneously via a multivari- 
ate GWAS approach. [47] This new strategy, which utilizes 
GWAS data more efficiently than univariate approaches, 
may help identify pleiotropic genes underlying diseases of 
shared genetic susceptibility, thereby revealing intercon- 
nected pathophysiological networks for a spectrum of 
common human diseases. 

The first bivariate GWAS analysis for osteoporosis and 
obesity, conducted in 1,000 unrelated US Caucasians, sug- 
gested the 50X6 gene's importance in co-regulation of os- 
teoporosis and obesity.[48] This gene was not significant 
at the genome-wide level in separate analyses for either 
BMD or BMI, highlighting the power of bivariate analyses 
in discovering novel pleotropic genes/loci. The SOX6 gene 
was independently found to be significant at the genome- 
wide level in GEFOS-2 [1 2] and was suggestively significant 
at the genome-wide level in the study of Hsu et al.[27] 
Identifying genes with pleiotropic effects, such as SOX6, 
may offer deeper insights into the pathophysiology of os- 
teoporosis, and potentially, may reveal key mechanistic 
links between osteoporosis and obesity, two major public 
health problems. 

5. The first whole genome sequencing study 

GWASs typically identify common variants important for 
common human disorders. Many susceptible rare variants 
may be missed in GWASs due to the limited power and/or 
resolution of genotyping.[49,50] Recent advances in next- 
generation sequencing technologies have greatly en- 
hanced our ability to discover functional rare variants.[51] 
Styrkarsdottir et al.[9] published the first whole genome 
sequencing study for osteoporosis. Through whole-ge- 



nome sequencing of 2,230 Icelandic individuals with dis- 
cordant BMD values, they found a rare nonsense mutation 
within the leucine-rich repeat containing G protein-cou- 
pled receptor 4 (LGR4) gene that was strongly associated 
with low BMD and OF. This mutation leads to termination 
of LGR4 at position 126 and fully disrupts its function. A 
potential problem with OF analysis is that all types of frac- 
tures were grouped together, which may introduce false 
positive or negative results as aforementioned. The pheno- 
type of carriers of the c.376C>T mutation overlaps that of 
Lgr4 mutant mice.[9] Interestingly, although this mutation 
was associated with a wide range of phenotypes across 
species (i.e., humans and mice), it was not present in Dan- 
ish and Australians.[9] Therefore, its effects in other human 
populations need to be further evaluated and validated. 

INTERPRETATION OF INCONSISTENT 
RESULTS 

Through the identification of novel genes, CNVs, and bi- 
ological pathways that are associated with osteoporosis, 
recent GWASs and large-scale meta-analyses have greatly 
advanced our understanding of the pathophysiology of 
osteoporosis and OF. As illustrated above, however, com- 
parison of the results across individual GWASs and meta- 
analyses raises a number of critical questions. 

1) Are individual GWASs still useful when meta-analysis, 
with a much larger sample size, is available? 

2) Why did some genes/loci identified as being signifi- 
cant in independent individual GWASs (e.g., ADAMTS18 in 
Ref [30]) fail to attain significance in independent GEFOS-2 
meta-analyses,[1 2] which had a much larger sample size? 

3) Why did some genes/loci identified as being signifi- 
cant in individual GWASs (e.g., RAP1A, TBC1D8, and OSBP- 
L1A in Ref [27]) fail to attain significance at the genome- 
wide level in GEFOS-2, in circumstances where the samples 
used in the individual GWASs were included in GEFOS-2 
meta-analyses?[12] 

4) Why have inconsistent findings even been observed 
between meta-analyses (e.g., GEFOS-1 [1 1 ] and GEFOS-2[1 2] 
) whose samples overlapped to a large extent? 

5) Should the meta-analysis with the largest number of 
individuals be considered as the gold standard to evaluate 
findings of other independent studies (especially relatively 
small GWASs)? 
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To address these critical questions, we performed a se- 
ries of theoretical analyses using simulation studies and 
published our results elsewhere.[52] Our theoretical analy- 
ses, under ideal situations, demonstrated that: 

1) Although the power of an individual GWAS study (of 
average sample size) to identify any particular locus (of av- 
erage effect size) is limited, the power to identify at least 
one (any one) locus can be high. This may explain the ob- 
servation that a number of previous individual GWASs 
have identified novel loci despite the relatively limited 
sample size of each of the studies. Given the anticipated 
large number of significant loci that have eluded detection 
thus far, individual studies are still valuable in identifying at 
least some of these underlying effect loci. 

2) The number of loci that can be detected in meta-anal- 
yses greatly exceeds the number that can be detected in 
individual GWASs. However, the power of a meta-analysis 
to identify many independent loci simultaneously can still 
be limited. 

3) The meta-analysis has rather limited power to repli- 
cate the findings of particular loci from individual GWASs 
at the genome-wide significance level, particularly for 
SNPs with small effects, implying inconsistent findings be- 
tween independent individual GWASs and meta-analysis is 
not unexpected. 

4) Adding heterogeneous samples into a subset of ho- 
mogeneous samples can reduce power for a meta-analysis, 
rather than having the anticipated effect of increasing 
power due to increased sample size. This was clearly shown 
in our published theoretical study[52] which examined the 
effects of sample heterogeneity by simulating a subset of 
samples with true effects and the other subset of samples 
without effects in a meta-analysis (based on GEFOS-2). Fig- 
ure (adapted from [52]) illustrates that heterogeneity re- 
sults in power loss in both the mixed-effects and random- 
effects models. This effect was further demonstrated in our 
more comprehensive theoretical analyses which considered 
a wide range of situations and scenarios.[53] Therefore, se- 
lecting samples with homogeneous effects may be at least 
as important as enlarging sample size for meta-analysis. 
This observation may help explain the inconsistent findings 
between GEFOS-1 [11] and GEFOS-2 [12] despite the fact 
that the latter included the majority of samples in the for- 
mer, and approximately 3 times as many total samples. 

Our theoretical analysis was not intended to denigrate 
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Fig. 1. Power of meta-analysis in heterogeneous populations. We 
simulated 17 studies (total sample size of 32,961 study subjects), 7 
studies having phenotypic effects and 10 studies having no pheno- 
typic effects. For simplicity, each sample was simulated with MAF of 
0.3. Between-study variance was set at 0.6. "Subset" samples were 
those having effects, and "Total" were total samples. "_F" and "_R" 
denote "fixed-effects" and "random-effects" models for meta-analy- 
sis. Significance level of meta-analysis was set at 5x1 0 s . Power was 
estimated based on 10,000 replications. 

the value and importance of GWAS meta-analysis. Large- 
scale meta-analysis clearly represents a powerful tool for 
identifying novel genetic variants for osteoporosis. Never- 
theless, it is critical to recognize that meta-analysis may in- 
troduce false positive/negative results due to heterogeneity 
and other confounding factors. Consequently, it is impor- 
tant to exert caution when designing studies, and interpret- 
ing results generated by meta-analysis of GWASs. In particu- 
lar, the results of meta-analysis should not be used to evalu- 
ate the validity of findings from individual GWASs. Well-de- 
signed individual GWASs, using homogenous samples, have 
significant potential to identify additional novel genes/loci 
and related biological pathways for osteoporosis and OF 
that may not necessarily be significant in meta-analysis with 
larger samples, of more variable heterogeneity. 

IMPLICATIONS AND PERSPECTIVES 

Over the past 4 years, individual GWASs and meta-analy- 
ses of GWASs have highlighted more than 70 genes/loci 
and related biological pathways that contribute to the 
pathophysiology of osteoporosis and/or OF. Some of these 
genes, such as WNT, RANK-RANKL-OPG, have known func- 
tional relevance to bone metabolism and endochondral 
ossification, and their contribution to osteoporosis has 
been well established in earlier candidate gene studies. 
[5,54,55] The function of these established genes, proteins 
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and related biological pathways has been reviewed else- 
where.[10] However, the purpose of GWASs is not merely 
to verify previous findings but, more importantly, to identi- 
fy novel genes and pathways associated with complex dis- 
eases. It is notable that more than half of the genes/loci 
listed in Table 1 are novel, though their functional impor- 
tance to bone metabolism awaits validation, ultimately 
through molecularfunctional studies. 

Collectively, the genes/loci identified from individual 
GWASs and meta-analyses, to date, explain less than 6% of 
the variance in BMD variation. Therefore, further endeavors 
are needed to explore undiscovered genetic factors associ- 
ated with BMD variation. At the DNA level, there are several 
possible paths towards uncovering these novel and elusive 
genes: 

1) Well-designed, powerful individual GWASs. For such 
studies to be powerful, well-defined and homogeneous 
phenotypes measured with high data quality and accuracy, 
large and homogeneous samples, and comprehensive sta- 
tistical and bioinformatical analyses should be needed. 

2) Meta-analyses with even larger samples than GE- 
FOS-2, which must be executed with caution because of 
the power loss associated with between-study heteroge- 
neity and other confounding factors. 

3) Meta-analyses with smaller sample sizes, but with less 
between-sample heterogeneity, to reduce the probability 
of generating false positive and false negative results. 

4) Utilizing available GWAS data to perform additional 
analyses such as CNV, pathway based, and multivariate 
analyses. 

5) Performing genetic studies focused directly on OF, 
rather than BMD or other less critical osteoporosis pheno- 
types. For OF studies, large and homogeneous samples of 
the same type of OF are critically important and essential. 
Although it is tempting to mix different types of OFs to in- 
crease sample size, heterogeneity may contribute to false 
findings and results that complicate data interpretation. 

Current evidence suggests that the genetic architecture 
of osteoporosis and OF is complex, involving both com- 
mon and rare functional variants.[50,51] These findings are 
similar for other complex human diseases that represent 
significant public health problems (e.g., obesity, diabetes, 
and cancers). GWASs and meta-analyses are largely de- 
signed to identify common variants. Thus, many suscepti- 
ble rare variants may be missed in GWASs due to the limit- 



ed power and/or resolution of genotyping. Recent advanc- 
es in next-generation sequencing technologies, however, 
have greatly enhanced our ability to discover functional 
rare variants.[51] Unfortunately, large-scale whole genome 
re-sequencing studies are currently prohibitively expen- 
sive. More economically feasible approaches might involve 
re-sequencing of targeted genes/loci, or exome sequenc- 
ing studies. Another potential approach to identify rare 
variants is to perform imputation analyses based on pub- 
licly available genome sequence data. The recent 1 ,000 Ge- 
nomes Project (http://www.1000genomes.org/) produced 
a comprehensive catalog of human genomic variants, in 
particular rare variants.[56,57] Thus, through imputation of 
currently available GWAS samples, it is now feasible to 
identify/infer the majority of known human genomic vari- 
ants (including rare variants) for these association studies. 
As technological advances decrease the unit price for ge- 
nome re-sequencing, it is expected that large and power- 
ful genome-wide re-sequencing studies will become feasi- 
ble in the near future, resulting in the identification of nu- 
merous rare genetic variants that effect osteoporosis and 
OF. It is critical to recognize that in order to identify rare 
variants with reasonable power, re-sequencing studies will 
require very large sample sizes due to the low minor allele 
frequencies of these rare variants. 

GWASs and re-sequencing studies focus on DNA se- 
quence variants alone. Recent studies, however, have 
shown that epigenetic regulation is also critical to the 
pathophysiology of many human complex diseases, and 
may partially account for the heritability that has not been 
accounted for based solely on DNA variants.[58] Epigenetics 
refers to reversible, heritable changes in gene regulation 
that occur without a change in DNA sequence.[59] Conse- 
quently, in order to identify additional and novel heritable 
factors contributing to BMD variation and osteoporosis risk, 
it has become necessary and timely to study epigenomic 
regulation in relevant bone-related cells, including DNA 
methylation, histone modification, microRNA and long non- 
coding RNAs. 

Protein post-translational modifications, such as protein 
phosphorylation, have also been shown to play significant 
roles in regulating gene expression,[60,61 ] signal transduc- 
tion, [62,63] and cellular functions[64] closely related to 
bone metabolism. Functional studies at the protein level, 
including protein expression[28,65] and post-translational 
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modification, may also contribute to our comprehensive 
identification and understanding of cellular and molecular 
mechanisms regulating bone metabolism. 

GWASs are ultimately studies at the DNA level, and the 
above comments illustrate the importance of studying 
complex human diseases from a systems biology perspec- 
tive.[66,67] Gene expression is a complex process that is 
regulated simultaneously and interactively at DNA, RNA, 
protein, epigenomic and environmental levels. Therefore, a 
genomic convergence or systems biology approach that 
organically integrates the information from GWASs, gene 
expression, proteomics, epigenomics, protein post-transla- 
tional modification, and gene-environment studies may 
help facilitate the identification of key pathways that are 
globally involved in the pathogenesis of osteoporosis and 
OF. For example, using a systems genetic analytic ap- 
proach, Calabrese et al. identified a physiologically relevant 
gene network and used it to discover novel genes and reg- 
ulatory mechanisms involved in the function of osteoblast- 
lineage cells.[68] In another study, Deng et al. ascertained 
SOD2 as a susceptibility gene for osteoporosis in Chinese 
by integrating evidence from DNA, RNA, and protein levels. 
[69] Ultimately, the functional relevance of the identified 
variants needs to be confirmed by in vivo and/or in vitro 
molecular biology studies. 
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